Chinese Journal of Astronomy and Astrophysics manuscript no. 
(Version: HuangYF.tex; printed on February 2, 2008; 1:16) 



Dynamical Evolution of Gamma-Ray Burst 
-<<y Remnants with Evolving Radiative Efficiency 

o : 

^ ! J. B. Feng 1 , Y. F. Huang 1 ' 2 , Z. G. Dai 1 and T. Lu 1 ' 2 



o 

o 



> 

in 
\o 
o 

(N 
O 

6 



1 Department of Astronomy, Nanjing University, Nanjing 210093, China 

2 LCRHEA, Institute for High-Energy Physics, Chinese Academy of Sciences, Beijing 
100039, China 

Received 2002 ; accepted 2002 

Abstract In previous works, a generic dynamical model has been suggested 
by Huang et al., which is shown to be correct for both adiabatic and radiative 
blastwaves, and in both ultra-relativistic and non-relativistic phases. In deriving 
their equations, Huang et al. have assumed that the radiative efficiency of the 
fireball is constant. They then applied their model directly to realistic cases where 
the radiative efficiency evolves with time. In this paper, we abandon the above 
assumption and re-derive a more accurate dynamical equation for gamma-ray 
burst remnants. Numerical results show that Huang et al.'s model is accurate 
enough in general cases. 
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1 INTRODUCTION 

Although the progenitors of gamma-ray bursts (GRBs) are still controversial (Cheng & 
Dai 2001; Cheng & Lu 2001b; Lu et al. 2000a, b), it is generally believed that energetic 
fireballs should be involved, where baryons are eventually accelerated to ultra-relativistic 
speed (Wu et al. 2001). After the main burst phase, the thin baryonic shell expands at 
ultra-relativistic speed into the surrounding matter, producing afterglows in soft bands 
(Cheng, Huang, & Lu 2001; Mao & Wang 2001a, b; Gou et al. 2001a, b; Huang, Yang, 
& Lu 2001; Zhang & Meszaros 2002). For good recent reviews on afterglow observations 
and theories, see van Paradijs, Kouveliotou, & Wijers (2000) and Cheng & Lu (2001a). 
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The dynamics of the gamma-ray burst remnants is different in two cases in which the 
remnant expansion is either adiabatic or highly radiative (Blandford & McKee 1976, 
1977). However, the conditions under which the remnant dynamics may be consid- 
ered adiabatic or radiative are far from unambiguous and are crucially dependent on 
poorly known questions about postshock energy exchange between protons and electrons 
(Meszaros, Rees, & Wijers 1998). Furthermore, a partially radiative regime with decreas- 
ing radiative efficiency may exist in realistic fireballs (Dai, Huang, & Lu 1999). So, it is 
necessary to construct a dynamical model that is able to describe a realistic fireball, i.e., 
a fireball with evolving radiative efficiency. 

The dynamics of gamma-ray burst remnants has been studied extensively (Sari 1997; 
Cohen, Piran, & Sari 1998; Panaitescu, Meszaros, & Rees 1998; Wei & Lu 1998; Chiang & 
Dermer 1999; Rhoads 1999; Panaitescu & Meszaros 1999; Kobayashi, Piran, & Sari 1999; 
Huang et al. 1998a, b, c, 1999a, b, 2000a, b, c, 2002; Dermer & Humi 2001). Especially, a 
generic dynamical model was proposed by Huang, Dai and Lu (1999a, hereafter HDL99), 
which is shown to be applicable to both ultra-relativistic and non-relativistic blastwaves, 
no matter whether they are adiabatic or highly radiative. In their derivation, Huang, Dai 
and Lu implicitly assumed that the radiative efficiency of the fireball, e, is a constant 
during the deceleration. They then generalized their model to discuss realistic blastwaves, 
where e evolves with time (Huang et al. 2000a, b). In this work we will inspect their 
generalization carefully. We first repeat the derivation of HDL99, but abandoning the 
e = const assumption. We then compare our result with that of HDL99 numerically. It is 
found that Huang et al.'s generic model can be applied to realistic remnants satisfactorily. 



We assume that after the initial GRB phase, the total energy left in the fireball is com- 
parable to the radiation energy emitted in gamma-rays, i.e., Eq ~ 10 51 -10 52 ergs. Denote 
the mass of the contaminating baryons as Mo, then the fireball continues to expand at 
a Lorentz factor of rj = Eq/(Mqc 2 ). Subsequently, at a radius Rq, the expansion of the 
fireball starts to be significantly influenced by the swept-up medium and external shock 
may form (Rees & Meszaros 1992). As usual, i?o is supposed to be 



where n is the number density of the interstellar medium, m p is the mass of a proton. 
2.1 Basic dynamical equations 

In HDL99, a generic dynamical model that is applicable in both ultra-relativistic and 
non-relativistic phases of GRB afterglows has been proposed. The key point of the model 



2 DYNAMICS 




(1) 



Dynamical Evolution of GRB Remnants with Evolving Radiative Efficiency 



3 



is a differential equation 




dm Mq + em + 2(1 — e)^m' 

where 7 is the bulk Lorentz factor of the fireball, m is the swept-up mass. Equation 
is derived as follows. Global conservation of energy implies that 



Here U is the co-moving internal energy with rest-mass excluded, J7 ra( j is the internal 
energy that is radiated from the fireball. If a fraction e of swept-up kinetic energy is 
instantaneously radiated from the fireball, then dU la d = — e(7 — \)dmc 2 . The internal 
energy U in the fireball changes because of the change of the kinetic energy of the swept- 
up matter, due to expansion of the fireball and the energy loss through radiation. Thus, 
we assume U = (1 — e)U cx , where U cx is the internal energy produced in this expansion. It 
is usually assumed that dU cx = (7 — l)dmc 2 . However, the jump conditions (Blandford & 
Mckee 1976) at the forward shock imply that U cx = (7— l)mc 2 , so the correct expression 
for dU cx under thin shell approximation should be dU cx — d[(j — l)mc 2 ] = (7 — l)dmc 2 + 
mc 2 d"f. Assuming e ~ const, then from equation (^) we can obtain equation (|J). 

It is worth noting that in the expression of dU cx = (7 — l)dmc 2 + mc 2 d r y, the term 
mc 2 dj is negative when the fireball is decelerating. This term, in fact, represents the loss 
of internal energy due to volume expansion of the fireball, i.e., the adiabatic loss term 
{dUadi) defined by Dermer and Humi (2001). This can be clearly seen from equation 
(13) of Dermer and Humi (2001). Under thin shell approximation, their equation can be 
approximately simplified as mc 2 d r y. 

In the above derivation, e is assumed to be constant during the deceleration. However, 
in realistic fireballs, e is expected to evolve from 1 to owing to the changes in the relative 
importance of synchrotron-induced and expansion-induced loss of energy (Dai, Huang, & 
Lu 1999). Equation (||) has been simply generalized to the case that e evolves with time 
(Huang et al. 1999b, 2000a, b, c). However, this might induce some errors. Below, we will 
abandon the constant e assumption and derive the equations that are strictly applicable 
for fireballs with evolving radiative efficiency. 

The assumption that U = (1 — e)U cx overestimates the true internal energy, because 
at late stages e is near 0, but at early stages it is about 1. Instead of using U = (1 — e)U ex , 
we use the expression that dU = (1 — e)dU cx . Substituting it into equation (|^), we obtain 
another differential equation depicting the evolution of fireballs 



d[7(M c 2 + mc 2 + U)] = dmc 2 + jdU lad . 



(3) 



d- t 




(4) 



dm 



M a + m + U/c 2 + (1 - e)jm ' 



with 



dU = (1 - e)dU ( 



ex 



(1 - e)[(7 - l)dmc 2 + mc 2 dj]. 



(5) 
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In the highly radiative case (e ~ 1, and U — 0), equation (Q) reduces to the case of 
Blandford & McKee (1976) 

*L = ^ Z 1 (6) 
dm M + m ' ' 

While in the fully adiabatic case (e ~ 0, and U = U ex = (7 — 1)toc 2 ), equation (^) reduces 

to the adiabatic case of HDL99 

*L = - . (7) 

dm Mo + 27m 

In fact, taking e = const, equation (||) exactly reduces to the generic model of HDL99. 
But if e evolves with time, we would expect that the fireball described by equation (^) 
will decelerate more rapidly than another fireball described by equation (||). 

2.2 Radiative Efficiency 

According to Blandford & McKee (1976), the electron number density (n!) and energy 
density (e') of the shocked medium in the frame co-moving with the fireball can be 
written as (also see: Huang et al. 1998b) 

, 77 + I 



7-1 

, 77 + 1 



(8) 



e 



— (7 - l)nm p c , (9) 



7- 

where 7 is the adiabatic index of the shocked medium, which is generally between 4/3 
and 5/3. Equations (||) and (|^) are appropriate for both relativistic and non-relativistic 
blastwaves. From the definition of 7 (Blandford & McKee 1976), Dai, Huang, & Lu (1999) 
gave a simple and useful approximate expression for 7: 7 ~ (47 + l)/(37). It can be seen 
from this approximation that 7 ~ 4/3 for an extremely relativistic blastwave and 7 ~ 5/3 
for a non-relativistic shock. 

As usual, we assume that the magnetic density in the co-moving frame is a fixed 
fraction €3 of the internal energy density, viz., B' = {STTtBG 1 ) 1 / 2 , and that the shock- 
accelerated electrons behind the blastwave carry a fraction e e of the internal energy 
(Huang et al. 2000a, b). This implies that the minimum Lorentz factor of the random 
motion of electrons in the co-moving frame is Jcmin = £e(7 — ^)m p /m e + 1. We here 
consider only synchrotron emission from these electrons, and neglect the contribution of 
inverse Compton emission because the latter emission is of minor importance particularly 
at late times of the evolution (Waxman 1997; Dai & Lu 1998). The energy of a typical 
accelerated electron behind the blastwave is lost both through synchrotron radiation and 
through expansion of the fireball, thus the radiative efficiency of this single electron is 
given by t^/(*^ + ^ x )(Dai & Lu 1998; Dai, Huang, & Lu 1999), where t' syn is the 
synchrotron cooling time, t' syn = 6irm e c/(o-TB' 2 "f eym i n ), and t' ex is the co-moving frame 
expansion time, t' ex — R/(-fc). Here R is the radius of the blastwave. Since all of the 
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Fig. 1 Evolution of the bulk Lorentz factor (7). The dashed line corresponds 
to Eq. (|^). The solid line is drawn according to Eq. (Q). Parameters: Eq = 10 52 
ergs, n = l cm" 3 , M = 2 x 1O _5 M , e e = 1.0 and e B = 0.01. 



accelerated electrons behind the blastwave carry only a fraction e e of the internal energy, 
the radiative efficiency of the fireball can be given by (Dai, Huang, & Lu 1999) 

At - (10) 



In the highly radiative case, e e ~ 1 and t' syn <C t' ex , we have e ~ 1. The early evolution of 
the remnants is likely in this regime. For an adiabatic expansion, e e <C 1 or t' syn 3> t' exl 
we get e ~ 0. The late evolution is believed to be in this regime. In realistic case, the 
radiative efficiency of the fireball (e) evolves from about 1 to (Huang et al. 2000a). 

2.3 Numerical Results 

The evolution of the radius and swept-up mass are described by (Huang et al. 1998a, 
2000a, b) 

dm = AirR 2 nm p dR, (11) 



dR = /3c 7 (7 + V7 2 - l)dt, (12) 

where t is the time measured in the observer's frame. Then equations (^|) and Q can be 
solved numerically. 

Figure 1 compares the evolution of the Lorentz factor calculated according to 
equations (0) and (0). In our calculations, we take Eq = 10 52 ergs, n = 1 cm~ 3 , 
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Fig. 2 Evolution of the shock radius (R). Parameters and line styles are the 
same as in Fig. 1. 




Fig. 3 Evolution of the radiative efficiency of the fireball (e). Parameters are 
the same as in Fig. 1. 
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Fig. 4 Evolution of the fireball momentum (P) under different assumptions 
for the radiative efficiency of the fireball (e) . The solid line is the case when e 
evolves according to equation (|l0|). The dashed line is the adiabatic case, i.e., 
e = 0. The dotted line is the highly radiative case, viz., e = 1. Parameters are 
the same as in Fig. 1. 



Mo = 2 x fCP 5 M Q , e e = 1.0, cb = 0.01. In both cases, equation ( JlOj ) is used to de- 
pict the evolution of e. We see that, as expected above, the bulk Lorentz factor of the 
fireball (7) calculated by equation (0) (the solid line) declines more rapidly than that of 
equation (the dashed line). But we notice that the difference is slight. Figure 2 shows 
the time dependence of the blastwave radius (R). Figure 3 shows the evolution of the 
radiative efficiency of the realistic fireball (e) . 

The relation between the radius (R) and the fireball momentum (P — — l) 1 ' 2 ) is 
shown in Figure 4. The solid line is the case when e evolves according to equation (|To|) . 
The dashed line is the adiabatic case, i.e., e = 0. The dotted line is the highly radiative 
case, viz., e = 1. We can see that, at early times when the realistic fireball is ultra- 
relativistic and highly radiative, the solid line approximately satisfies P oc R~ 3 . At late 
times when the fireball is non-relativistic and adiabatic, the deceleration is approximately 
P oc i?~ 3 / 2 , consistent with the Sedov limit. 

We emphasize that for the e = const cases, the results are precisely the same in the 
two models characterized by equation ^ and equation (|]). 
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Fig. 5 Predicted afterglow light curves in fixed frequency v = 10 15 Hz. S v is 
in units of ergs s _1 cm -2 Hz -1 . The dashed line corresponds to the generic 
model of HDL99, and the solid line corresponds to Eq. ([I]). Parameters adopted: 
E Q = 10 52 ergs, n = 1 cm~ 3 , M = 2xl(T 5 A/ Q , e e = 1.0, e B = 0.01, p = 2.1, and 
D = 1 Gpc. Note that after ~ 10 4 s, slopes of the two curves are approximately 
identical. 



3 LIGHT CURVE 

In section 2, the dynamical evolution of a postburst fireball has been calculated numeri- 
cally. As in Dai Huang & Lu (1999), we calculate the light curves of optical afterglows. 
The results are shown in Figure 5. Here the solid line is drawn by using the dynamics 
of equation (0) and the dashed line is drawn by using equation (||). We see that the 
difference between the two curves is not notable. Flux densities on the dashed curve are 
higher by about 2 after the light curve peak, but the slopes of the two curves are identi- 
cal. Please note that in our calculation, we have taken relatively large parameter values 
for e e and es- £e = 1.0, eg = 0.01. If these two parameters are taken typical values as 
e e ~ 0.1, 6b ~ 10 -4 — 10 -6 , then the difference will be even smaller. 

4 DISCUSSION AND CONCLUSIONS 

The generic model of HDL99 is applicable to both radiative and adiabatic fireballs, and 
in both ultra-relativistic and non-relativistic phases. A possible problem is that whether 
this model is correct or not when the radiative efficiency of the blastwave (e) evolves 
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with time. We have shown that in this case, for the evolution of 7 and R, the errors 
induced by the generic model is almost negligible. Errors in the optical light curves are 
slightly amplified due to the strong dependence of flux density on the Lorentz factor, but 
the results are still acceptable. We suggest that the generic model in its simple form of 
equation (7) in HDL99 could be safely used when e varies during the deceleration. 

A dynamical model that is applicable to both relativistic and non-relativistic expan- 
sion has been established for quasars and active galactic nuclei by Blandford & McKee 
(1977). Their dynamics is most convenient for either adiabatic or highly radiative blast- 
waves, even allowing for steady injection of energy into the remnant from the central 
engine. However, for partially radiative blastwaves, especially blastwaves with an evolv- 
ing efficiency, the simple generic dynamical model of HDL99 is still more convenient. 
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